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Abstract 



We investigate the out of equilibrium dynamics of global chiral supersym- 
metry at finite energy density. We concentrate on two specific models. The 
first is the massive Wess-Zumino model which we study in a selfconsistent 
one-loop approximation. We find that for energy densities above a certain 
threshold, the fields are driven dynamically to a point in field space at which 
the fermionic component of the superfield is massless. The state, however is 
found to be unstable, indicating a breakdown of the one-loop approximation. 
To investigate further, we consider an O(N) massive chiral model which is 
solved exactly in the large N limit. For sufficiently high energy densities, we 
find that for late times the fields reach a nonperturbative minimum of the ef- 
fective potential degenerate with the perturbative minimum. This minimum 
is a true attractor for 0(N) invariant states at high energy densities, and 
this provides a mechanism for determining which of the otherwise degenerate 
vacua is chosen by the dynamics. The final state for large energy density is a 
cloud of massless particles (both bosons and fermions) around this new non- 
perturbative supersymmetric minimum. By introducing boson masses which 
softly break the supersymmetry, we demonstrate a see-saw mechanism for 
generating small fermion masses. We discuss some of the cosmological impli- 
cations of our results. 

PACS numbers: ll.10.-z, 11.30.Pb, 12.60.Jv, ll.15.Pg 
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I. INTRODUCTION 



Supersymmetry M, by nature of its Grassman transformation parameter, behaves differ- 
ently from ordinary symmetries under the effects of finite temperature, as was first recognized 
over 20 years ago by Das and Kaku || and has been studied since by a number of authors 
f3HT3|. In particular, it was shown that unbroken supersymmetry in equilibrium at zero 
temperature becomes broken at finite temperatures 0J§]. As usual, such breaking of the 
continuous symmetry is necessarily signaled by the appearance of a massless particle, the 
Goldstone. Only in this case the Goldstone is a fermion rather than a boson, due to the 
fact that the symmetry transformation parameter is a Grassman variable; the residual flat 
direction in field space must correspond to a fermionic degree of freedom. The Goldstone 
fermion is referred to as the Goldstino. 

The fundamental reason why supersymmetry becomes broken at finite temperature is 
easy to understand. At finite temperature, fermions and bosons obey different statistics. 
This means that while the Lagrangian may admit a supersymmetry between boson and 
fermion components, the way in which these components are populated at finite temperature 
breaks any such symmetry. In equilibrium, the only state which transforms supersymmetri- 
cally is the zero temperature supersymmetric ground state. 

These issues are relevant for a number of reasons. We see no supersymmetric partners to 
the fermions nor to the bosons in nature. Hence, if we are to assume that supersymmetry is 
fundamental, then it must be broken. Supersymmetry breaking is an appealing possibility 
which raises the question of the nature and identification of the resulting Goldstino. 

The possible breaking of supersymmetry is also relevant to the early universe. In particu- 
lar, many popular inflation models are based on supersymmetry in some form [14]. However, 
inflation occurs very far from thermal equilibrium and while supersymmetric models are un- 
derstood at zero and finite temperature, the non-equilibrium dynamics of such models have 
yet to be properly studied. And even when supersymmetric models are considered in the 
context of inflation, it is common practice to discard some of the degrees of freedom as 
irrelevant. However, if the results in thermal equilibrium are to be a guide, this may be a 
dangerous thing to do. Processes such as the breaking of supersymmetry and the consequent 
appearance of massless degrees of freedom might also occur out of equilibrium as the avail- 
able energy is distributed differently among fermions and bosons. This may be particularly 
relevant for the process of reheating after inflation, for which light fermions may play an 



important role [15-19 



We examine these issues by means of explicit numerical solutions of supersymmetric toy 
models allowed to evolve far from equilibrium. First results describing the dynamics of O(N) 
chiral supersymmetry at finite energy were published recently [M. Although the primary 
applications of interest are cosmological, we simplify the analysis by restricting ourselves to 
Minkowski spacetime. 

We begin our analysis with a study of the massive Wess-Zumino model in a self-consistent 
one-loop approximation. We find two distinct dynamical regimes. At sufficiently low energy 
density, the system remains near one of the two supersymmetric vacua of the model. The 
non-zero energy is taken up in the oscillation of the field zero modes about their perturbative 
vacuum values. Some particles are produced, but there is no obvious massless fermion in 
the spectrum over the lifetime of the numerical simulations. These states are continuously 
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connected to the zero temperature vacuum states in the sense that as one reduces the initial 
energy density of the system toward zero, the system approaches the zero temperature 
vacuum without any changes in the qualitative behavior of the system. 

The second regime is quite different. When the system exceeds a certain critical energy 
density such that the fields can sufficiently sample the region separating the two perturbative 
minima, the system becomes driven to the point halfway between the two minima, i.e., to 
the local maximum of the effective potential. What is special about this point is that the 
mass of the fermion field precisely vanishes. 

There are a couple things we learn from this behavior. First, this provides another 
instance of the dangers of relying solely on equilibrium constructs such as the effective 
potential when dealing with systems far from equilibrium The perturbative effective 

potential alone would never cause one to expect that the system would approach such a 
point. Second, with sufficient energy to reach the relevant point in field space, the system 
will invariably approach a state for which the fermion field is massless. Numerical tests of 
the system indicate that this behavior is obtained for arbitrarily high energy densities. 

However, we also find that this state has an instability to this order in perturbation 
theory as it includes a scalar field with negative mass squared, corresponding to a so-called 
spinodal instability. This is expected from being at the local maximum of the effective 
potential. This is a signal of the breakdown of perturbation theory and an indication that 
the state which eventually forms is non-perturbative in nature. This unfortunately limits 
our study in the one component model to the early time behavior before the instability 
becomes important. 

In order to proceed beyond the perturbative Wess-Zumino model, we introduce a model 
with an internal O(N) symmetry, which can be solved exactly, even into the non-perturbative 
regime, in the limit iV — > oo. While an O(oo) field theory may not be realistic, it provides a 
useful self-consistent toy model for the study of possible characteristics of realistic theories 
- particularly those with continuous symmetries - and may also be considered as a first 
approximation to systems with moderate values of N. Such techniques have a long history 
in field theory [^3 in general as well as in supersymmetric theories In the case of purely 



scalar models, large N studies provide concrete examples of the non-perturbative symmetry 
breaking processes leading to Goldstone bosons |p2| , ^5| and to a dynamical formation of a 
flat potential related to the Maxwell construction of equilibrium thermodynamics |26|] . 

Our results again indicate two regimes for states invariant under the O(N) symmetry. 
As with the one-loop Wess-Zumino model, there is a continuous spectrum of states at suffi- 
ciently low energy densities which are connected to the perturbative O(N) symmetric vacua. 
For sufficiently high energy densities, the system again is driven to a point in field space for 
which the fermion modes are massless. This point, <p = —m/X which acts as an attractor 
to the dynamics is a nonperturbative minimum of the effective potential degenerate with 
the perturbative minimum. The final state obtained by real time evolution for large en- 
ergy density is a cloud of massless particles (both bosons and fermions) around this new 
nonperturbative minimum. 

There are several interesting facts about this system that deserve mention. First, we 
see that the dynamics effectively chooses the vacuum for the system during the high energy 
density phase. In a cosmological context, once such a vacuum is chosen, the universe would 
stay in that vacuum indefinitely. We therefore see a mechanism for choosing one vacuum 
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over another in spite of the fact that they are degenerate in energy. Second, we find that all 
particles are massless for arbitrarily high energy density. This is to be contrasted with the 
general, non-supersymmetric case for which continuous symmetries are always restored at 
high enough energy density and particles become massive. The study of the spontaneously 
broken O(N) scalar theory provides a good example of this more usual behavior [p6| . In the 
present case, however, the supersymmetry acts to protect the kind of Goldstone phase we 
find with all particles, bosons and fermions, being massless. 

As a final example, we examine what occurs when the supersymmetry of the O(N) 
model is softly broken by the introduction of additional small scalar mass terms. The result 
is that while there is still a set of massless bosons in the high energy density phase, their 
superpartner bosons and fermions gain masses. While the bosons gain a mass proportional 
to the soft breaking mass scale, there is a see-saw mechanism for the fermions which gives 
them a mass proportional to the square of the soft breaking mass scale divided by the 
overall scale of supersymmetry, providing a very natural mechanism for producing very light 
fermions which could be relevant to neutrino mass generation. 

We continue in the next section with the introduction of the Wess-Zumino model and 
provide the renormalized one-loop equations of motion appropriate to an out of equilibrium 
study. This is followed by our numerical results for the model, showing the two distinct 
regimes at low and high energies. In Section [TTT| , we introduce the O(N) model and again 
provide our numerical results. In Section |iy], we include soft supersymmetry breaking into 
the O(N) Lagrangian and show how the various fields gain non-zero masses. Our conclusions 
are provided in Section [V]. We also present two appendices with details of the renormalization 
procedures used in the non-equilibrium formalism. 



II. THE WESS-ZUMINO MODEL TO ONE LOOP ORDER 
A. model and equations of motion 



We consider the supersymmetric Wess-Zumino model [[27]]. It is based on a single chiral 
super-multiplet S with superpotential JTJ 

W(S) = ^mS ■ S + \\S ■ S ■ S , (2.1) 

This can be broken down into component fields via S = (A, B; if?; F, G), where A is a scalar, 
B is a pseudo-scalar, ip is a Majorana fermion, and F and G are scalar and pseudo-scalar 
auxiliary fields respectively. After eliminating the auxiliary fields, the Lagrangian is given 

by 



C = \d,Ad»A + \d,Bd»B - \m 2 (A 2 + B 2 ) - \\mA (A 2 + B 2 ) - h 2 (A 2 + B 2 
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2 p 2 

- -mijip - -Aipip - % —Bi)^i\) . (2.2) 

In order to follow the dynamics, it is convenient to break up the scalar field A into its 
expectation value and small fluctuations about that value: 
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A = 4> + A , <p= (A) , (A) = . (2.3) 

For convenience, we take B to have zero expectation value. We write 

B = B , (B) = 0. (2.4) 

Likewise, we treat ip as a pure fluctuation with (ip) = 0. We can then expand the Lagrangian 
in orders of the fluctuations, A, B, and ip. 
In zeroth order we find 

£<W = - V{4>) , (2.5) 

with the classical potential 

V(<f>) = \^ 2 (m + ^ , (2.6) 

as sketched in Fig. 1. 




FIG. 1: The classical scalar potential with A = m = 1.0. The dotted line represents 
the spinodal line above which there is sufficient energy to enter the negative curvature 
portion of the potential, a = — m(l + l/y/S)/X, b = — m(l — l/\/3)/A. The upper dashed 
line indicates the initial value for which the energy is high enough for cp to pass the 
maximum, c = m 4 /8A, the lower dashed line marks the upper limit for the initial value 
for (j), which leads to a stable configuration, d = m 4 /18A. Initial values in the region 
above d lead to unstable configurations. 

This yields the classical part of the equation of motion. The first order in the fluctu- 
ation vanishes because the expectation value of the fluctuations is zero. The second order 
expression reads 

£(2) = l - d ,Ad»A + \d»Bd»B ~^(A 2 + B 2 ) 

~m\<P (3A 2 + B 2 ) - |V (3A 2 + B 2 ) 

i 1 A - 

+-ip0i/> - -mipip - -<jnpip . (2.7) 
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We can derive the equation of motion for the classical field and for the fluctuations from 
these Lagrangians. We find for the classical field in one-loop approximation: 



+ V'{<j>) + ^ (m + A0) [3(A 2 ) + (B 2 )} + ~$V> = 



while for the fluctuations we find 



A [ - V 2 + m 2 + -A 2 2 + 3mA0 j A = , 

D ( -V 2 + m 2 + ^A 2 2 + mA^) 5 = 0, 
(i^ — m — \<fi) ij) — . 



(2.8) 

(2.9) 

(2.10) 
(2.11) 



It is convenient to collect the mass term and the terms depending on <p(t) into time- 
dependent masses 



M 2 A (t) =m 2 + 3mA0(t) + -A 2 2 (t) 

A 2 

M%(t) =m 2 + m\(f)(t) + — 2 (t) , 



M^{t) =m + X<f>(t) . 
We note that the supersymmetry relation 



Ml+M 



2M; 



o 



(2.12) 

(2.13) 
(2.14) 



(2.15) 



is satisfied by the time-dependent masses for all times. We further note that in the region 

(2.16) 



-l--L<A^/™<-l + -L, 



M. 2 A is negative, leading to instabilities. M. 2 B (t) on the other hand is everywhere positive 
definite. 

We provide integral expressions for the expectation values appearing in Eq. ( p.8|) below. 
Note that if the fermion mass Ai^(t) vanishes, i.e. <p(t) is equal — m/A, the contribution 
of the bosonic fields A and B cancels in ( |2.8|) . In this case, only the fermionic fluctuations 
influence the equation of motion and the behavior of the classical field. 

We expand the scalar fields in terms of mode functions f A and and the fermion field 
in terms of the spinor solutions of the Dirac equation as in [To*| . |2"S| 



A{t,x) 
B{t,x) 
V(t,x) 



d 3 k 


1 


(2?r) 3 


2ujao 


d 3 k 


1 


(2tt) 3 


2ubo 


d 3 k 


1 


(2tt) 3 


2uj^q 



CA, Mt) + C AM f* A (t) 

CB, Mt)+C^ k fUt) 

b k , s U Ks (t)+bK s V- k>s (t) 



ikx 



ikx 



(2.17) 
(2.18) 
(2.19) 



with the usual (anti-)commutation relations for the time independent annihilation and cre- 
ation operators 
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[c jM , 4 >k ,] = 2^ (27T) 3 5 3 (k - k') , j — A,B 
{KsA>, s } = 2^o(2vr) 3 5 3 (k- k')5 ss , . 

The frequencies in Eqs. 17| ) are defined as 
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k 2 + m 2 j0 ,j = A,B,i/> 



(2.20) 
(2.21) 

(2.22) 



with rrij = _Mj(0). 

The equations of motion for the bosonic mode functions are given by 

f j (t)+[k 2 + M 2 j (t)]f j (t)=0, 

with j = A, B. The initial conditions for the fields are chosen as 

/,(0) = 1, ./•(()) 



(2.23) 



(2.24) 



For the fermions, we define mode functions f^{k,t) and g^(k,t) through the relations 



U s (k,t) = N [id t + H k (t)) U(k,t) 



Xs 

, o 



K(k, t) = N [id t + W_ k (t)] g^k, t) 



'0 

\Xs 



(2.25) 
(2.26) 



where the \s with s = ±1 are helicity eigenstates with eigenvalue s and TL^it) is defined as 

?4(t) = 7o7 • k + • (2.27) 

The mode functions f^{k,t) obey the second order differential equation 



^ - iM^t) + k 2 + Ml(t) 



U(k,t) = 0, 



while g^(k,t) = f^(k,t). The initial conditions are 

U(k, 0) = 1, f^k, 0) = -iu^o . 

The integrals appearing in the equation of motion for 0, Eq. ( |2.8| ), are: 

d 3 k \f A {t)\ 2 
(2tt)3 2u A0 ' 
d 3 k \f B (t)\ 2 



(A 2 ){t) - 
(B 2 )(t) - 

<W(*) = - 



(2tt) 3 2lu B0 
d 3 k 1 



(27r) 3 CJ^o 



2uj 



2k 2 



l/*(*)f 



(2.28) 
(2.29) 

(2.30) 
(2.31) 
(2.32) 



The energy density can be calculated as the trace over the Hamiltonian. With the results 
for the scalar fields in [^8] and the fermion fields in |H| we can express the energy density 
in terms of the mode functions in the following way 
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£ = \i + VW)] 



+ 
+ 
+ 



d 3 k 


1 


(2tt) 3 




d 3 k 


1 


(2tt) 3 




d 3 k 


1 



(2tt) 3 2^ 



~\\fA\ 2 + \{k* + M\{t))\f A \- 
'l\f B \ 2 + l(k* + Ml(t))\f B \ 



(2.33) 



It is easy to verify that the time derivative of the energy vanishes if the equation of motion 
for and those of the fluctuation modes are satisfied. 

Both the equation of motion for the classical field and for the energy density are divergent. 
Therefore we have to consider their renormalization. Within a computational scheme based 
and extended in |28[] for non-equilibrium dynamics it is possible to make a clean 
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separation between divergent and finite parts, allowing us to directly compute the finite 
parts of the integrals appearing in eqs.( |2.30| )-( |2.32| ). Details of the scheme, developed for 
various physical models in refs. |T5|j28|j30| -p5| , are presented in the Appendix A. 

The structure of the renormalized equations is found to be analogous to that of the 
unrenormalized ones, so we do not display them here. The divergent expectation values 
appearing in Eqs. (|2.30|) - (|2.33|) are replaced by finite, subtracted expressions. The resulting 
divergent counter terms are the same as the standard ones obtained in equilibrium. It is 
worth mentioning that supersymmetry leads to the expected cancellations of quadratic and 
quartic divergences leaving only the relatively well behaved logarithmic divergences - the 
primary advantage of supersymmetry models. 



B. Results and interpretation 

We find two distinct phases possible in the dynamics of the Wess-Zumino model. We 
examine each in turn. 





FIG. 2: Zero mode evolution showing the 
low energy density phase for m = A = 1.0 
and 0(0) = 0.23. 



FIG. 3: The effective squared masses; 
solid line: Ml(t), dotted line: M 2 A (t), 
dashed line: M. 2 B {t)\ the parameters as in 
Fig. 2. 



S 



At low energy density the model behaves much as one would expect from the scalar 
potential shown in Fig. 1. The zero mode <fi oscillates about one of the two vacua, see 
Fig. 2. Particle production is minimal. As a result of the deviation of the zero mode from 
the vacuum, the masses of the various quanta each oscillate about their supersymmetric 
values as shown in Fig. 3. Note that while there is a small mass splitting between the 
masses, the overall supersymmetry sum rule, ( [2.1 5|) remains satisfied. As the overall energy 
density is reduced to zero, these mass splittings vanish, indicating that one reaches the true 
supersymmetric vacuum. 

These numerical results indicate that the system goes to a limiting cycle. Namely, the 
zero mode keeps oscillating forever. This is in contrast to the $ 4 model in either the large 
N limit or in the Hartree approximation for which the zero mode always has a constant 
infinite time limit [[^,[36[ . 

The behavior of the system is very different, however, if the initial energy density is 
large enough. In particular, we find that whenever the initial energy density is sufficiently 
large, there is an attractor solution which causes the system to fall into a state signaled 
dramatically by the vanishing of the fermion mass. Examination of the classical potential 
for <p (Fig. 1) and the expression for the effective mass of the fermion %p ( |2.14| ) reveals that 
the only point for which the fermion is massless corresponds precisely to the local maximum 
of the classical potential for 0! 




10 20 30 

t 



FIG. 4: Zero mode evolution showing the 
high energy density phase for m = 1.0, 
A = 2.0 and 0(0) = 0.5. 
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FIG. 5: The effective squared masses; 
solid line: M^(t), dotted line: M\(t), 
dashed line: M. 2 B (t); parameters as in 
Fig. 4. Note that -> while M\ 

becomes negative. 



The results of a sample evolution is plotted in Fig. 4 with the effective masses of the 
fields given in Fig. 5. We find, indeed, that 4> goes to the point — m/A with the mass of 
the fermion approaching zero. We also see that the evolution starts to blow up toward the 
end of the simulations as demonstrated in Fig. 6 (and this unstable behavior continues if 
allowed to evolve further). This is due simply to the fact that <ft is sitting on top of its 
potential, leading to a severe instability in its fluctuations (A 2 ) due to the effective negative 
mass squared. This is also the signal of the breakdown of the one loop approximation. In 
fact, it is expected that perturbation theory as a whole will break down at this point due to 
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the process known as spinodal decomposition |57j], the result being the formation of a fully 
non-perturbative and inhomogeneous state. 
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-0.5 - 
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-1.1 — 
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Total Energy 



FIG. 6: The scalar fluctuation (A 2 ) show- 
ing unstable behavior, parameters as in 
Fig. 4. 



FIG. 7: Asymptotic value of as a func- 
tion of the total energy with m = A = 
1.0. Notice, that 0^ = — m/A = —1.0 
since the energy density is larger than the 
threshold. 



Although the dynamics breaks down at one loop order, we can nevertheless consider what 
will happen as the evolution proceeds into the non-perturbative regime. We are helped by 
the sum rule ( |2.15| ) for the masses of the fields which states that A4 2 A + M. 2 B = IM. 2 ^ (we 
mention once again that this relation holds both for out of equilibrium unbroken and broken 
super symmetry). First, we note that the one loop results satisfy this relation, only with the 
caveat that the squared mass for A is negative. Experience with non-perturbative dynamics 
Pr]| , |57| -pU| in purely scalar field theories provides clues as to what to expect. The main 
point is that the resulting growth of fluctuations of A due to the instability will provide a 
contribution tending to increase the effective mass of the field toward zero (and at the same 
time decreasing the B field mass toward zero). The end result will be massless fields A and 
B. This is, in fact, the only possible stable solution of the sum rule given massless fermions. 

We therefore expect that the concave portion of the classical potential will become flat- 
tened out by non-perturbative effects due to a spinodal instability in the field A. This 
flattening has been seen in studies of scalar A$ 4 models in the large A" limit of an O(N) 
field |38j and using the self-consistent Hartree approximation for a single scalar field |39| . 



These studies verified that the flattening of the potential is related to the dynamical approach 
to the Maxwell construction of the true free energy of the system [p6| , [4(i|1 . 

It is worthwhile to discuss the processes that lead to the transition between the two 
dynamical regimes. 

If the initial energy is sufficient for the system entering the spinodal region, Eq. ( |2 . 1 6|) , 
the region between a and b displayed in Fig. 1, the instability of the low momentum modes 
of the field A comes into play. This spinodal energy density is given by p s = m 4 /18A 2 , 
marked by d in Fig. 1, and if we start the motion with 0(0) > 0, the corresponding initial 



amplitude is S = m(— 1 + y5/3)/A. If 0(0) exceeds S only slightly the system enters the 
unstable region only for a short time, and the instability does not build up. At late times 
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the classical field has transferred energy to the quantum modes and no longer enters the 
spinodal region; it again ends up in the regime of stationary oscillation. 

At higher excitations the time spent by the system in the spinodal region increases and 
after a few oscillations the system is trapped there. The unstable low momentum modes of 
the field A grow indefinitely and the classical amplitude tends to 0oo = —m/X, i. e., at the 
maximum of the classical potential. In the phase diagram in Fig. 7 we have displayed the 
final value of in dependence of the total energy. The first part of the curve indicates that 
for low energies 0^ ends up near the minimum of the classical potential as described above. 
At some initial value 0(0) = cr ; t we find a rather sharp transition into the unstable regime, 
0oo = -m/X. 

This dynamics continues for energies above the energy density of the maximum of the 
potential Vm ax = m 4 /8A, and to the highest excitations we have considered. Above p = Vm ax 
the oscillations of the system extend into the region of the second supersymmetric minimum 
= —2m/ A and back again to = 0. With each oscillation the system enters the spinodal 
region and the instability can build up uninhibited. The result is that the system again ends 
up at 0oo = —m/X. This final state is not a stationary state in the true sense, however; 
while the classical amplitude approaches = —m/X the fluctuations of the field A grow 
indefinitely, see Fig. 6. 



C. Fields near the local maximum of the potential 



We can analytically solve the field equations fl2.8|) , fl2.17|) and (|2.23|) for the order pa- 
rameter near the local maximum of the potential : 

We set, 



m 

' x ■ 



0(t) = -- + A(i) 

and we will restrict to times where A(t) <C m. 
Eqs.( |2.23"D for thus becomes 



(2.34) 



hit) + 



2 2 w 



f A {t) = 



-h 2 



showing that the A-modes with k < m/\[2 are growing exponentially as e V 2 * . There- 
fore, the quantum fluctuations (A 2 ) [see eq.( |2.30D l will be dominated by these spinodally 
unstable low- A; modes and will grow as 



(2tt) 3 2<j A0 



(m 



t)3/2" 



where the function D(t) is the order one for late times. The equation of motion for the zero 
mode (12.81) takes the form 



A(t) + 1-^ + 1^ [A 2 W + 3(A 2 )(t) + (B 2 )(t)} | A(t) + \ X^m) = • (2.35) 
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Since as t grows the zero mode approaches — y, A(t) — > and (A 2 )(t) — > oo. Therefore, in 
order eq. (|2.35|) keep valid, the product 



(A 2 )(t) A(t) = G(t) 



must stay bounded. Moreover, it must fulfill 

1 



G(t) 



3A 



[2.36) 



a relation which is satisfied by the dynamics in the numerical simulations, see Fig. 8. 
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FIG. 8: Numerical check of (A 2 ) A = -(^)/3A, dashed line: (A 2 ) A, dotted line: 
-W)/3A. 



III. LARGE N EXTENSION OF THE WESS-ZUMINO MODEL 

Having seen that the self-consistent one-loop approximation to the Wess-Zumino model 
breaks down due to the spinodal instability, we now turn to an extension of the model 
for which we can follow the evolution into the non-perturbative domain. We couple the 
Wess-Zumino superfield to iV chiral superfields satisfying an internal O(N) symmetry. By 
taking the large N limit, we arrive at a semi-classical supersymmetric model for which we 
can numerically solve for the exact and fully non-perturbative field dynamics. 

A. Model and equations of motion 



Our O(N) extension of the Wess-Zumino model consists of a chiral superfield multiplet 

50 = (A , B ; ipo] Fq, Go), which acts as a singlet under 0(N), coupled to iV chiral superfields 

51 = (Ai, Bf, ipi] Fi, Gj) with i = 1 . . . N and which transform as a vector under O(N). The 
superpotential has the form 



C = -MSo ■ So 



K 



1 N A 
zSn • So ■ So -|- —771 / Sj ■ Sj -\- — 1= 7 So ■ Si • Sj . 



N 



8=1 



(3.1) 



i=i 
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Notice that the model considered in Ref. |41] is a special case of ours for m = k = 0. 
Again, we expand in terms of the component fields and eliminate the auxiliary fields via 
their equations of motion. To allow for a consistent large N limit, the expectation value of 
A is taken to be of order y/N. We set 

(A ) = v/iV0 , (B Q ) = . (3.2) 

The latter condition, which amounts to a choice of initial conditions, is chosen because it 
significantly simplifies the equations of motion. We assume that the initial state satisfies 
the O(N) symmetry which requires that (Aj) = (Bi) = 0. Given the O(N) symmetry, it 
is convenient to define fields A, B, and ip such that J2iAiAi = NA 2 , J2iB{Bi = NB 2 , and 
J^i^i^i — Nipip. In taking the large N limit, particular care must be given to terms of the 
form J2i AiBi which turn out to be of order \/N . This is most easily seen by computing 
the squared quantity J2iJ2j AiBiAjBj] the resulting delta function 5j contributes a factor 
1/N which cancels one of the factors of iV coming from the summations. The Lagrangian 
to leading order in N is 



£ _ i„ 

N ~ 2 

1 

2 



+ -daAd^A + -d.Bd^B - -<p 2 (m+ 
2 2 M 2 V 2 



— (m 



+ \<pf {A 2 

1 



B' 



A 



M0 + -k^ 2 + -A (A 2 - B 2 )] (A 2 - B 



-mipip X(j)ifjifj . 



(3.3) 



Note that the fluctuations of the supersymmetry multiplet So do not appear to this order 
in N, with only the mean field <fi contributing to the Lagrangian. The consequence is that 
the Lagrangian of this theory is equivalent to a supersymmetric O(N) field 5^ coupled to a 
fully classical background zero mode <fi. This is a consequence of the semi-classical nature 
of the large iV approximation. 

The equation of motion for the mean field is found to be 



+ (M + K(, 



M(j)+ \ K ^ + \X(A 2 -B 2 ) 



+ A(m + \<P)(A 2 + B 2 ) + -A(^) = , 

while the fermion obeys 

[i$ - (m + \<j>)] ip = , 
and the mode functions from which the scalar fluctuations are built satisfy 

U + k 2 f A +(m + \<P) 2 f A + A [M0 + 1 -k<\> 2 + l -\(A 2 - B 2 )j f A = 
h + k 2 f B + (m + X<P) 2 f B - A [M0 + \k$ 2 + \\(A 2 - B 2 )] f B = 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



The integrals for the fluctuations are constructed just as in the ordinary Wess-Zumino model 
via Eqs. fl2.30|) -( p.32j ). It is convenient to define the masses 
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M 2 _ = X 
such that 



M^ = m + \(p, (3.8) 

(3.9) 



M0+^ 2 + \\(A 2 -B 2 ) 



M 2 A = Ml + M 2 _, (3.10) 
M 2 B = M%- M 2 _ . (3.11) 

We see again that the supersymmetry sum rule A4 A + M. 2 B — 2A4 2 = is automatically 
satisfied. 

The equation for M 2 , (|3.9| ), plays the role of a gap equation expected in the large N 
framework, since the fluctuations appearing in the integrals on the right hand side depend 
upon the masses A4 A and Ai b which in turn depend upon M._. This becomes more explicit 
in the fully renormalized form found in Appendix B. We also note that M. 2 _ can be either 
positive or negative, with the consequence that either or both A4 A and A4 2 B can become 
negative during the evolution. 

The energy density is 

£ = ^ + i(^ + ^ 

+i(i 2 + k 2 A 2 + M\{t)A 2 ) + ^(B 2 + k 2 B 2 + M 2 B (t)B 2 } (3.12) 

k2 



Again the equations of motion and the energy density have to be renormalized. Despite 
the non-perturbative nature of the large N limit, one of its important properties is that it 
is possible to consistently renormalize the theory. The details are presented in Appendix B. 
We separate out the divergent terms so that they may be treated analytically, leaving the 
finite parts to be included in numerical simulations. This allows us the freedom to choose a 
regularization scheme without regard to constraints of the numerical simulations. We then 
use dimensional regularization to define the counterterms. Again the structure of the finite 
equations remains essentially the same as the one of the bare equations, with divergent 
expectation values replaced by finite ones, and with finite corrections to the masses and 
couplings. We find consistency with the renormalization of the large- iV equilibrium theory. 



B. The supersymmetry in the large N limit 

Examining the leading order Lagrangian, Eq. ( |3.3|) , the question arises as to whether 
taking the large N limit is consistent with supersymmetry. One might be particularly 
concerned since the singlet superfield Sq appearing in ( |3.1|) is represented only by the single 
scalar zero mode = (Aq) / yN and one might expect that it would be necessary to have a 
corresponding fermion field into which may transform. To clarify this point we examine 
the transformation properties of Aq. Under supersymmetry, Aq transforms as 

4, = A + 5A = A o + C0o , (3.13) 
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where £ is the Grassman supersymmetry transformation parameter. Taking the expectation 
value of Eq. (|3.13| ) yields the transformation law for 0: 

0' = (A' )/VN = (Ao)/VN + (C^o)/ViV = + C(^o)/ViV = , (3.14) 

where the last equality follows from the vanishing expectation value of ipo ■ We therefore 
see that is invariant under supersymmetry transformations. This is consistent with the 
treatment of as a classical background zero mode coupled to the supersymmetric O(N) 
multiplet. 

To leading order, it is no longer necessary to consider the transformation properties of 
the O(N) singlet fields corresponding to the superfield So as none of these fields appear in 
our large N Lagrangian fl3.3| ). The remaining transformations are: 



50 = 0, (3.15) 

5A = (y, (3.16) 

5B = ClBip, (3.17) 

5i[> = (i0 + m + A0) (A + 7s S) ( , (3. 18) 

6$ = C [-i0 (A - l5 B) + (m + A0) (A + l5 B)] . (3.19) 



Through use of the equations of motion for A, B, and ■0, it is straightforward to show that 
the variation of the Lagrangian ( |3.3| ) vanishes under these supersymmetry transformations 
up to a total derivative. Therefore, we see that to this order, the Lagrangian is completely 
supersymmetric with acting as a classical background field. 

C. Results and interpretation 

As in the ordinary Wess-Zumino model, we find two distinct phases, one at low energy 
densities and one at high energy densities. 

Unlike the ordinary Wess-Zumino model, however, this model has perturbative vacua for 
which the fermion masses vanish. 

In the high energy density phase, the time evolution leads to a final state formed by a 
cloud of massless particles (both bosons and fermions) around a minimum = — m/X of the 
effective potential, degenerate with the tree level minimum. This nonperturbative minimum 
defines a new sector of the theory. [The perturbative sector is formed by particles around 
the perturbative ground state = 0.] Since this nonperturbative minimum at = —m/X 
is degenerate with the perturbative minimum, it is invariant under supersymmetry. Hence, 
supersymmetry is not broken at this nonperturbative minimum. In addition, the fact that 
all masses vanish at the new minimum = —m/X supports the supersymmetric character 
of this point. In addition, we choose O(N) invariant states since {A/) = (Bj) = for all 
times. 

However, the final state obtained by real time evolution for large energy density is formed 
by a cloud of massless particles (both bosons and fermions) around this new nonperturba- 
tive minimum = —m/X. This state indeed breaks supersymmetry since the energy is 
distributed differently among the fermions and bosons on the top of this zero-energy non- 
perturbative ground state due to their differing statistics. This is analogous to which is 
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known from equilibrium studies at finite temperature where a thermal gas of particles is 
around a perturbative supersymmetric vacuum. 

In summary, for large energy density the system goes to a nonperturbative minimum 
at = —m/A which acts as an attractor to the dynamics and where a kind of Goldstone 
phase develops with all particles, bosons and fermions, being massless. Interestingly, these 
phenomena happens for the highest energy densities tested. The reason that keeps all 
particles massless is that supersymmetry results in cancellations in the contributions to the 
effective masses of the fields, see eqs. (|3.9| ) - (|3.11| ). 



The situation here looks similar to the case where there is no symmetry restoration at 



high energies and massless Goldstone bosons are therefore present P2| , fi3|| . While individual 
field fluctuations may become large as the energy is increased, the net contribution of the 
fluctuations does not grow. Thus, arbitrarily high energy densities need not yield the mass 
corrections which ordinarily lead to symmetry restoration in spontaneously broken O(N)- 
theories pfi] . 

We now examine each phase in detail. 

The low energy, massive phase occurs in much the same way as in the ordinary Wess- 
Zumino model and is depicted in Fig. 9. with the corresponding masses shown in Fig. 10. 
However, because of the more complicated structure of the potential in the large N model, 
the range of energy densities for which the low energy density phase persists depend on the 
initial conditions. One can distinguish two cases: 

1. 0o > —M/k and m/A < M/k, or O < —M/k and m/A > M/k. In either of 
these two cases, the zero mode begins on the same side of its potential as the new 
nonperturbative supersymmetric minimum. This results in a relatively low critical 
energy density for reaching the high energy density phase. 

2. 0o > —M/k, and m/A > M/k, or 0o < —M/k and m/A < M/k. In these cases, 
there is an additional potential barrier between the initial value of and the new 
nonperturbative supersymmetric minimum. As a result, the transition to the high 
energy density phase occurs at a higher critical value of the energy density. 




50 100 150 200 

t 



FIG. 9: Zero mode evolution showing the 
low energy density phase for m = M = 
1.0, k = A = 1.0 and 0(0) = 0.2. 




50 100 150 200 



FIG. 10: The effective squared masses; 
solid line: Ml(t), dotted line: M 2 A (t), 
dashed line: M%(t); the parameters as in 
Fig. 9. 
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1 

FIG. 11: Zero mode evolution showing the 
low energy density phase for m = 1.0, 
M = 4.0, k = A = 1.0 and 0(0) = 0.2. 



20 40 60 80 100 

t 

FIG. 12: The effective squared masses; 
solid line: Ml(t), dotted line: M\{t), 
dashed line: M. 2 B (t); parameters as in 
Fig. 11. 



A second example of low energy density evolution is shown in Figs. 11 and 12 for which 
the masses are such that may efficiently decay into A and B particles, leading to apparent 
dissipation. We see that 0> decays and settles at a point near the classical minimum at = 0. 
The shift of the finite density minimum from this vacuum value is due to the growth of the 
fluctuations of the fields A, B, and ip, which are also responsible for the deviation of the 
masses from the supersymmetry value m. 

An example of the high energy density phase is depicted in Fig. 13 with the masses in 
Fig. 14. The evolution begins with large oscillations of over the entire classically allowed 
range of evolution. During this initial period, the field fluctuations (A 2 ) and (B 2 ) grow. 
After a relatively short period of time, the mean field settles down precisely to the point 
= —m/X. The result is that the N fermions become massless. As in the one-loop Wess 
Zumino case, this is an attractor state at high energy densities. However, in the present 
model the massless state is stable. A look to the effective potential |20| shows that 0> = — m/X 
is a zero-energy minimum of the effective potential. This minimum is degenerate with the 
perturbative vacuum = 0. 




FIG. 13: Zero mode evolution showing the 
high energy density phase for m = 1.0, 
M = 4.0, k = 2.0, A = 3.0 and 0(0) = 
-1.5. 



40 



30 -'i 




FIG. 14: The effective squared masses; 
solid line: -M^(t), dotted line: M 2 A (t), 
dashed line: Ai 2 B (t); parameters as in 
Fig. 13. 
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We therefore find that the system reaches a non-perturbative state for which each of 
the O(N) fermion and boson fields is massless and for which the vacuum energy identically 
vanishes as well indicating the presence of a new non-perturbative supersymmetric minimum 
at = —m/X. The fact that the massless state is an attractor is of great importance if one 
considers the model from a cosmological perspective. What it indicates is that the initial 
evolution of the system in the early universe can determine the ultimate vacuum state of 
the system, providing an effective means of selection between vacua which are otherwise 
degenerate in energy. 



IV. ADDING SOFT BREAKING TERMS 



One question that might be raised is what are the consequences of adding additional mass 
terms to the model which softly break the supersymmetry. This is particularly interesting 
because, as we have seen, the dynamics leads to massless fermions and, in the case of the 
O(N) model, massless bosons as well. Through soft supersymmetry breaking, it may be 
possible to introduce small masses to these final state particles, which would be determined 
not by the scale of supersymmetry, rather by some lower scale (e.g.., the electroweak scale) 
at which the soft supersymmetry breaking terms arise. 

To provide an example case, we introduce soft supersymmetry breaking to the O(N) 
model via a scalar mass m s for the A and B fields such that Eqs. ( 3.10|) and (|3.11 ) become 



M 2 A = M'^ + Mi + m 2 and M 2 B = M\ — M 2 _ + rrij. Such terms break the supersymmetry 
explicitly, while not producing any dangerous (i.e. non-logarithmic) divergences. We plot 
the results for the case m s /m = 0.1 in Fig. 15. 




FIG. 15: The effective field masses in the high energy density phase including soft 
breaking masses for the fields A and B; solid line: M.\(t), dotted line: M\(t), dashed 
line: M B (t); the parameters are m = 1.0, M = 4.0, m s = 0.1, k = 2.0, A = 3.0, 
(f)(0) = 0.9. 

The first thing to note is that we retain the high energy density attractor solution. As 
the new terms do not explicitly break the O(N) symmetry, there is necessarily one set of 
asymptotically massless Goldstone bosons, in this case represented by the B field. However, 
these fields' superpartners are no longer massless as they pick up contributions proportional 
to the soft supersymmetry breaking mass scale m s . For A, this mass is given approximately 
by A4a — v^rrisi corresponding to A4_ ~ m s and A4^ ~ 0. For the fermions, there is a see- 
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saw mechanism producing the mass such that the leading order contribution is proportional 
to m 2 s /m. This produces fermion masses which are suppressed by a factor of m s /m relative 
to their massive bosonic superpartners. 

V. DISCUSSION AND CONCLUSIONS 

Let us first summarize our main results. 

1. At low energy density, both the one loop and the large N Wess-Zumino models are 
found to allow for finite density non-equilibrium states based on supersymmetric vacua. 
Finite density corrections break the mass degeneracy of the scalars and fermions some- 
what, but do not lead to obvious massless fermions. In the large N model, such states 
are stable, albeit time dependent. They are also stable in the ordinary Wess-Zumino 
model to one loop order. 

2. At high energy density, the dynamics of the Wess-Zumino model and its large N 
extension leads to massless fermions. These states are stable in the large N model, 
but are highly unstable at one loop order in the ordinary Wess-Zumino model. 

3. In the large N model, the high energy density phase is characterized for arbitrarily 
high energy densities by all of the O(N) vector fields, scalars and fermions, becoming 
massless, and by the vanishing of the total vacuum energy. The phase is formed by a 
cloud of massless particles (both bosons and fermions) around a new nonperturbative 
supersymmetric minimum at = —m/X. 

4. The introduction of explicit soft supersymmetry breaking terms to the large N La- 
grangian results in masses for some of the O(N) vector fields. In particular, we find 
that a soft supersymmetry mass of m s for the scalars induces a mass for the A field 
M.A = V2m s and a fermion mass of order m 2 s /m, while B remains massless. 

We can visualize a number of consequences for cosmology and particle physics. 

1. It appears that we should expect massless fermions to occur in supersymmetric models 
of the early universe independent of whether the situation is one of equilibrium or far 
from equilibrium, or at least, if the universe is far from equilibrium, we should expect 
it to approach such a state. 

2. If the vacuum energy is indeed zero as in the large N case, then we have a convenient 
mechanism by which otherwise massive fermions (and possibly other fields as well) 
become massless. This could have the effect of priming the system so that soft super- 
symmetry breaking terms can yield fields with masses much smaller than the overall 
scale of supersymmetry. 

3. A final consequence of potential interest is that supersymmetry may protect continuous 
symmetries from being restored in the very early universe, as in our O(N) model. Such 
a result might alleviate the monopole problem, as without a symmetry restored phase, 
there would be no phase transition to produce such objects. 



19 



We have only begun to examine these possibilities, but many of these ideas can be pursued 
further using current techniques in out of equilibrium quantum field theory. One avenue of 
approach is to examine the ordinary Wess-Zumino beyond one loop order by means of mean 
field theory. This may clarify what happens to the masses of the A and B fields as well as 
having the potential to determine the vacuum energy of the resulting non-perturbative state. 
Also of interest is the behavior of these models in an expanding universe, as it would be 
important to see if a phase transition from the high energy phase to the low energy density 
phase occurs when energy is drained from the system via expansion. The relation between 
the expansion time scale and the relaxation time scale in the high energy density phase may 
also play a significant role. 

Of course, it would be beneficial to examine the detailed behavior of gauge multiplets 
in order to begin to get a better understanding of more realistic models of particle physics. 
It would also be very interesting to examine the possibility that neutrino masses could 
be realistically generated via an analogous mechanism to that of fermion mass generation 
in the softly broken supersymmetry model studied here. The techniques are available to 
tackle such problems, and the results of the present study encourage us to believe that 
there is much more to be learned from further studies of the out of equilibrium dynamics of 
super symmetric particle physics. 
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APPENDIX A: RENORMALIZATION: THE WESS ZUMINO MODEL 



In order to renormalize the one-loop-equations we use the perturbation scheme of 28 



which allows us to extract the divergences from the leading orders. We recall some basic 
equations and refer the interested reader to [0,0 f° r further details. We rewrite the mode 
equations (|2.23|) and ( p.28| ) as 



d 2 

dF 2 ~ 



with j = A,B,ip and 



(Al) 



V A (t) 


= 1" 2 (^ 2 


- <Pl) + 3mA (0 - 


<t>o) = M 2 A (t) - m 2 A0 


v B (t) 




- 0q) + m ^ (4> - 


j> ) = M 2 B (t) - m|o . 




= A 2 (0 2 - 


(pi) + 2\m (0 - <i 


i ) - «A0 




= M%{t)- 







(A2) 
(A3) 

(A4) 

With the initial conditions ( |2.24| ) it is possible to write the differential equation ( |A~T|) as an 
equivalent integral equation 
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— fdt' sin [u j0 (t - t')} V 3 {t')f 3 {t') • 
OOjQ Jo 



We can now make a general ansatz for the mode functions via 



fi(t) 



[1 + hM . 



(AS) 



(A6) 



Inserting this ansatz into the integral equation one derives an iteration that results in an 
expansion of the functions hj(t) in orders of the potentials Vj(t). 

The perturbative expansion for the bosonic fluctuation integrals leads to [EH1 



(A 2 ) = /-1K40) - V A (t)I- 3 (m A o) + (A\ Q , 
(B 2 ) = J_i(m B0 ) - V B (t)I- 3 (mBo) + (5 2 > fin • 

The integrals k — 1,3 are divergent. In dimensional regularization they become 

d 3 k 1 1 
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'2 , Anfi 2 
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d 3 £; 1 
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The finite parts are defined as 

d 3 k 1 f 
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(B 2 ) 



fin 



(2tt) 3 2^o 
d 3 fc 1 



\f A (t)\ 2 
\f B (t)\ 2 



2uj 2 m 



167T 2 ' 



v A (t) 



(A7) 
(A8) 



(A9) 
(A10) 



1 + 



2uj 2 m 



V B {t) 



(All) 
(A12) 



The integration over the subtracted integrand is finite. Note that using the ansatz (|A6|) 
and working with the functions hj(t) the leading divergence cancels explicitly, i.e., it is not 
included from the outset. 

For the fermions we obtain 



(#) = /- 3 (m 2 ) [2A^ + 4M 3 



T J t 



fin > 



with 



7&® 



16tt 2 



Am^M^t) - In 
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00 
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fin • 



(A13) 



(A14) 



It should be mentioned that the finite parts as defined here generally become singular as 
t — ► 0. These initial singularities can be removed by a Bogoliubov transformation ||44|| . We 
do not discuss this here, as it is not essential in the present context. 

The total contribution of the fluctuations in the equation of motion fl2.8|) has the form 



-\ 2 <P 3(A 2 ) + (B 2 ) 



+ 



mX 



3(A 2 ) + (B 2 ) + -(^) 
m 



(A15) 



Inserting the expressions for the expectation values given above we find 
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T = X 2 
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(A16) 



where J-fi n is defined by ( |A15| ), with all the expectation values replaced by their finite 
parts. Note that instead of the divergent integrals J_fc(mjo) with the different initial masses 
one has been able to collect the divergent part into one integral J_ 3 (m) where m is the 
renormalization mass, chosen here as the common physical mass parameter of the model. 
The dependence on the initial conditions appears in the finite terms 
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(A17) 
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The divergent part of (|A16|) can be written in the form 



F i n f A 



It can be removed by adding to 

1 



0-V'(0)]/_ 3 (m) . 

) a counter term superlagrangian 
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The equivalence between the mass and the coupling constant counter term is a special 
feature of supersymmetry. With these counter terms taken into account the equation of 
motion becomes finite. Explicitly it is given by 



1 + AZ) + (m 2 + Am)0 + -mA(l + AA)0 2 
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In the same way we have split T into divergent and finite parts we can decompose the 
fluctuation part of the energy as 

£&(t) = £fl,div(*) + £fl,fln(t) > ( A27 ) 

with 
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The energy has no initial singularity In addition to the quadratic and logarithmic diver- 
gences we find here a quartic one. In dimensional regularization it can be written as 
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In the same way as for the equation of motion we can now evaluate the momentum integrals 
in £a\ v , fix the counter terms, and find the finite contributions. This leads to 
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No divergent "cosmological constant" counter term is needed, as to be expected in a super- 
symmetric model. After taking into account the counter term Lagrangian ( |A22| ) with the 
previously determined coefficients the energy becomes finite and reads 



i(l + AZ)0 2 + l -{m 2 + Am)0 2 + -mA(l + AA)0 3 
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APPENDIX B: RENORMALIZATION OF THE LARGE- N MODEL 



The expressions for the expectation values are analogous to those of the simple Wess- 
Zumino model, and so is their perturbative expansion. Of course now the fluctuation masses 
A4j(t), given by Eqs. (|3.10|) , (|3.11|) and (|3.9|) , depend on the expectation values of the 
fluctuations, and so do the potentials Vj(t), defined in analogy to (|A2|)-([A~4|). 

In the following we will need the sum and the difference of the bosonic fluctuations. For 
the sum we obtain 



with 



1 



16vr 2 



In 



m 



AO 



M\{t) + \n 



m 



BO 



M 2 B (t)- m 2 A0 - m 2 m 



+ (A 2 ) Rn + (B 2 ) {in . 

(A 2 ) fi n and (-B 2 )fi n are again defined by ( All ) and ( |A12j ). For the difference we find 

1 . 1 



(A 2 - B 2 ) = -2/_ 3 (m 2 ) 



MX(f) + |kA0 2 + -X 2 (A 2 - B 2 ) 



with 



(Bl) 



(B2) 



(B3) 
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M\(t) 



m 12a 



M%{t) 



AO 



+ m 



BO 



+(A 2 ) fin - (B\ n . 



(B4) 



Obviously the equation for (A 2 — B 2 ) is implicit, as to be expected in the large- N framework. 
The logarithms of the ratios m 2 /m 2 arise from replacing the masses rrij by the common 
renormalized mass m of all the component fields in the integral I-^{m 2 ). It is convenient to 



rewrite these equations in terms of M._, as introduced in Eq. 



{A 2 - B 2 ) 



I -Am 2 



16tt 2 



In 



( m A0 m B0 \ 
\ m? ) 



. We get 



(B5) 



with 



171 AO ~ m B0 



+ (A 2 ) fin -(S 2 ) 



fill 



(B6) 



The expectation value of the fermionic fluctuations again decomposes as ( |A13| ) and ( |A14j ). 
We now introduce multiplicative renormalization factors by replacing 



A 

A - 



>Z X \ 
Z A A 



(B7) 
(B8) 
(B9) 
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and similarly for all other quantities. We first note that there is no divergent term that would 
require an infinite renormalization of the mass m. This mass occurs in the combination 
m + \<fi which now gets replaced by m + ZxZ^Xcp. We conclude that also the second term 
stays unrenormalized, and so 



We next analyze the mass M. 2 _{t). We have 

Ml (t) = Z M Z x Z^MX(f> + ^Z K Z x ZlK\<j) 2 + ^Z 2 A Z 2 X X 2 (A 2 - B 2 ) 
We have assumed Z& = Zb- Further, from (|B5|), we have 



m: 



l + \ 2 Z 2 A Z 2 x I„ 3 (m 2 ] 



Z M Z x Z^MX(f> + 1 -Z k Z x ZIk\$ 2 



(BIO) 



(Bll) 



(B12) 



+X 2 Z 2 A Z 2 



A n X 



In 



fmAom B o\ 



16tt 2 V 



m 



2 J 



M 2 _ + -Z 2 A Z\\ 2 T^ 



In order to get analogous relations for the finite quantities as for the unrenormalized ones 
we require 



1 + \ 2 Z 2 A Z\l^{m 2 ) 



z\z\ = z 



M 



Z k Za 



having used Eq. QBIOP . The first relation yields 

1 



rv2 rjl 

A X 



1 - A 2 /_ 3 (m 2 ) 



We then have further 



Z k Za 



1 



1 - A 2 /_ 3 (m 2 ) ' 
The renormalized equation for the mass M.- becomes 

A 2 



(B13) 



(B14) 



(B15) 



M 2 _(t) 



In 



fm A0 m B o\ 



or 



with 



16tt 2 V 



M 2 _{t) = C 



C 



MA0(t) + ^A0 2 (t) + \\ 2 T^ , 



M\ct>{t) + l -K\^{t) + l -\ 2 ^ n 



1 



A^ 



■In 



/ m A0 m B0 \ 
167T 2 \ m? J 

Using Eq. ([B10|) we also have the renormalized relations 



(B16) 



(B17) 



(B18) 
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M\{t) = [m + \<f){t)] 2 + M 2 _(t) 
M 2 B (t) = [m + A0(t)f -M 2 _(t) 
M^(t) = m + X<p(t) . 



(B19) 
(B20) 
(B21) 



We now turn to the equation of motion for the mean field Eq. ( |3.4|) . Introducing 
the renormalization factors we obtain 



+ -^-(M + K<j))M 2 _ + ZlZ^Zl^I^m 2 )^ 
-Z X Z 2 A X (m + A0) (A 2 + B 2 ) + Z X ZZ-(*V) = 



(B22) 
(B23) 



We first consider the last two terms. These terms are quadratically divergent. This "tadpole" 
contribution necessitates a cancellation between bosons and fermions. This works only if we 
postulate 



Zxf, — Za — Zb 



a relation indeed required by supersymmetry. Then 
Z X Z 2 A X [m + \<j>) (A 2 + B 2 ) + Z x Z 2 ^m) = Z X Z 2 A X 



M^A 2 + B 2 ) + \{^) 



(B24) 



(B25) 



Inserting the decomposition into divergent parts and fluctuation integrals the expression in 
brackets yields 



so that the equation of motion for <j) takes the form 

Z^ + |^(M + k<P)M 2 _ + Z{Z^Z 2 A \ 2 I^(m 2 
+ l -Z x Z\\FL + Z x Z 2 A X(m + \</>)f+ = . 

The coefficients of the two terms containing <ft combine into 

i + zizix 2 !.^ 2 : 1 - 1 



T - X 2 I,(m 2 ) 



(B26) 

(B27) 
(B28) 

(B29) 



Here we have used Eq. (|B14j ). The factor of the term agrees with that of the second term 
in Eq. QB27D 



Zm 



ZmZ$ 

Likewise for the fourth term in Eq. 



1 



l-X 2 lJm 2 ) 



we obtain 



X Zj a — Zj^Zj^Zj a — Zi^ 



1 - X 2 I^(m 2 ) 



(B30) 



(B31) 
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The same factor is obtained for the fifth term, so the renormalized equation reads 
+ \{M + kcI>)M 2 _ + + A(m + \<f>)?+ = . 



(B32) 



The various factors Zj are not yet completely determined. This is to be expected, since 
the equations of motion do not depend on the absolute normalization of the Lagrangian. 
When including the Z factors and using the various relations derived so far, the energy 
density is given by 



£ 



lzU 2 + lz 2 M Zl(M<p + 



+^Z 2 A (A 2 + k 2 A 2 + M 2 A (t)A 2 ) + ±Z 2 A (B 2 + k 2 B 2 + M 2 B {t)B 2 ) 



(B33) 
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The bosonic fluctuation energies defined clS clS 

£ A = ^(A 2 + k 2 A 2 + M 2 A (t)A 2 ) 
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(2tt)3 2uo A0 1 2 

and analogously for the field B. Its expansion in terms of divergent integrals reads 
£a = \h{m 2 M ) + h^(m 2 A0 )V A (t) - jI- 3 (m 2 A0 )V 2 A (t) + £ AM 



(B35) 
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Here we have used Eqs. (|A9|) , ( |A1U| ), and ( |A30[ ). The finite part is defined by subtracting 
the quantity 
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AAiv 
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(B36) 



from the expression ( |B34| ). 

The fermionic fluctuation energy is defined as 



£<< 



2 

d 3 k 1 



(B37) 



(2tt) 3 2^o 
It is expanded as 



[i [u^o - m^o] (f(k,t)f*(k,t) - f(k,t)f*(k,t)) - 2uj^M4t)} 
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(B38) 



where S^^ n is defined by subtracting 
d 3 k 1 



S, 



i/>,div 



(2tt) 3 2^o 
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(Ml 



m 



4ca 



</>0 2 



(B39) 



from the integral of Eq. (|B3j 

With these preparations we are ready to discuss the divergences of the energy density. 

The term proportional to M\ = \ 2 <p 2 adds to the kinetic term Z^(p 2 . The total factor 
(up to a finite term, see below) is Z 2 + Z\\ 2 I^. If we compare to the relations for the Zj 
found in the previous section we find that the prefactor of <j) 2 becomes finite if we choose 
Za = 1 and as a consequence Z\ — 1 — A 2 /_3. 

The second and the fifth term in the general expression for the energy, Eq. ( |B33|) can 
be rewritten, using Eq. ( Pll|) as 
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(B40) 



We note that ZmZ^/Zx = 1 for the choice ( |B43| ) , so that the first term on the right hand 
side is finite. The second term —M 4 i/2ZlX 2 has a factor Z^ 2 = 1 — A 2 /_3. This combines 
with the infinite parts of Sa, &b and £^ as 



--M 4 
2 - 



-/- 3 K 



--(Ml + M 2 ) 7_ 3 (m 
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(B41) 



So all infinite terms have cancelled, and the choice of renormalization constants 
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(B42) 
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(B44) 



based on the analysis of the equations of motion and of the energy momentum tensor. It 
remains to collect the finite terms. We find 
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(B45) 



with 



28 



£am — £a,&ti + 



1 



1287T 2 

1 

1287T 2 

1 

64^ 



m A0 - 4m 2 A0 M 2 A + 2 In 



m 



in 



Bo-4m| 7W| + 21n 



m 



A0 ^M A 



in 



Am 2 M Ml + 2\n 



m 



m 



-f )M% 

m z J 



(B46) 
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We note again that no 'cosmological constant' counter term has to be introduced. 
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